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Exotic dark matter together with the vacuum energy (associated with the cos- 
mological constant) seem to dominate in the Universe. Thus its direct detection 
is central to particle physics and cosmology. Supersymmetry provides a natural 
dark matter candidate, the lightest supersymmetric particle (LSP). Furthermore 
from the knowledge of the density and velocity distribution of the LSP, the quark 
substructure of the nucleon and the nuclear structure (form factor and/or spin 
response function), one is able to evaluate the event rate for LSP-nuclcus elastic 
scattering. The thus obtained event rates are, however, very low. So it is impera- 
tive to exploit the modulation effect, i.e. the dependence of the event rate on the 
Earth's motion and the directional signature, i.e. the dependence of the rate on 
the direction of the recoiling nucleus. In this paper we study such experimental 
signatures employing a supersymmetric model with universal boundary conditions 
at large tan/3. 
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I. Introduction 

In recent years the consideration of exotic dark matter has become necessary 
in order to close the Universe i The COBE data B suggest that CDM (Cold 
Dark Matter) is at least 60% a. On the other-hand evidence from two different 
teams, t^a-JIigh-z Supernova Search Team □ and the Supernova Cosmology 
Project □ '19 suggests that the Universe may be dominated by the cosmological 
constant A. Thus the situation can be adequately described by a baryonic 
component tts = 0.1 along with the exotic components flcDM — 0.3 and 
A = 0.6 (see next section for the definitions). In another analysis Turner 
gives J7 m = £Icdm + ^* rf^ 0-4. Since the non exotic component cannot 
exceed 40% of the CDM 0,0, there is room for the exotic WIMP's (Weakly 
Interacting Massive Particles). In fact the DAMA experiment □ has claimed 
the observation of one signal in direct detection of a WIMP, which wjtt 
statistics has subsequently been interpreted as a modulation signal 

In the most favored scenario of supersymmetry the LSP can be simply de- 
scribed as a Majorana fermionj-aJ|i|H**«i combination of the neutral components 
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of the gauginos and HiggsinostHl 



II. An Overview of Direct Detection - The Allowed SUSY Pa- 
rameter Space. 



Since this particle is expected to be very massive, m x > 30GeV, and extremely 
non relativistjc with average kinetic energy T < lOOKeV, it can be directly 
detected E3'E£l mainly via the recoiling of a nucleus (A,Z) in elastic scattering. 
In order to compute the event rate one needs the following ingredients: 

1) An effective Lagrangian at the elementary particle (quark) lej«eL ob- 
tained in the framework of supersymmetry as described , e.g., in Refs.EEJ. 

2) A procedure in going from the quark to the nucleon level, i.e. a quark 
model for the nucleon. The results depend crucially on the content of the 
nucleon in quarks other than u and d. This is particularly true for the scalar 
couplings as well as the isoscalar axial coupling l 18 l r°l , 

3) Compute the relevant nuclear matrix elements E3Ej using as reliable as 
possible many body nuclear wave functions. The situation is a bit simpler in 
the case of the scalar coupling, in which case one only needs the nuclear form 
factor. 

Since the obtained rates are very low, one would like to be able to exploit 
lulation of the event rates due to the earth's revolution around the sun 
this end one adopts a folding procedure assuming some distribu- 
I of velocities for the LSP. One also would like to know the directional 
rates, by observing the nucleus in a certain direction, which correlate with the 
tion of the sun around the center of the galaxy and the motion of the Earth 
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The calculation of this cross section has become pretty standard. One 
starts with representative inpu t in the restricted SUSY parameter space as 
described in the literature li3ii-3. We will adopt a phenomelogical procedure 
taking universal soft SUSY breaking terms at Mqjjt, i.e., a common mass for 
all scalar fields too, a common gaugino mass Myn and a common trilinear 
scalar coupling Aq , which we put equal to zero (we will discuss later the influ- 
ence of nonjjzero Ao's). Our effective theory below Mqut then depends on the 
parameters E3: 



m , Mi/2, Mo, "G, Mgut, h, , h, , h T , tan/3 , 



where o>q = 9qI^ (<?g being the GUT gauge coupling constant) and ht, hb, h T 
are respectively the top, bottom and tau Yukawa coupling constants at Mqjjt- 
The values of qq and Mqjjt are obtained as described in ReltU. For a specified 
value of tan f3 at Ms, we determine ht at Mqjjt by fixing the top quark mass 
at the center of its experimental range, mtim t ) = 166GeV. The value of h T 
at Mqjjt is fixed by using the running tau lepton mass at mz, m T (mz) — 
1.746GeV. The value of hb at Mqjjt used is such that: 
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m b {m z )sM = 2.90 ± 0.14 GeV. 

after including the SUSY threshold correction. The SUSY parameter space is 
subject to the following constraints: 

1. ) The LSP relic abundance will satisfy the cosmological constrain: 

0.09 < fl LSP h 2 < 0.22 (1) 

2. ) The Higgs bound obtained from recent CDF@ and LEP20, i.e. m h > 113 GeV. 

3. ) We will limit ourselves to LSP-nucleon cross sections for the scalar cou- 
pling, which gives detectable rates 

4 x io- 7 P b < < 2 x 10- 5 P b (2) 

We should remember that the event rate does not depend only on the nucleon 
cross section, but on other parameters also, mainly on the LSP mass and the 
nucleus used in target. The condition on the nucleon cross section imposes 
severe constraints on the acceptable parameter space. In particular in our 
model it restricts tanfi to values tanfi ~ p5fl. We will not elaborate further on 
this point, since it has already appeared til. 



III. Expressions for the Differential Cross Section . 



The effectivfi-Lagrangian describing the LSP-nucleus cross section can be cast 
in the formEij 

Leff = -^|{(xi7 A 75Xi)^A + (xiXi)J} (3) 

where 

Jx = N^(f°+f^r 3 + f Al5 + f\ l5 T 3 )N , J — N (/ 6 ° + f}r 3 )N (4) 

We have neglected the uninteresting pseudoscalar and tensor currents. 
Note that, due to the Majorana nature of the LSP, Xi7 A Xi — (identically). 

je above ingredients the differential cross section can be cast in the 

form 
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<M "-'' ] = 2{^bvf [{ts + tv F2 ( u )+^ pin F n (u)] (.H 



- Vr{A) ,2 r A 2 r/ f f 1 A-2Z 2 g 2 fl r (A) 2 
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spin j- 

-'spin — ^ n v o Sspin i sspin f o u[U) If I 

3(1 + 4) 2 

J A 



s(u) = [Sn o{ o)f^fl + 2 4no(o)nx(o)^ + ^i(o)) 2 ] (8) 



= < x o Cv (9) 
A - 2Z^ 2u v 0s2 1 2r)+l (2u 



(1 + 4) 2 fv A 



(2^6)2 (l + ry) 2 (v*) 1 
(10) 



cr* x0 = proton cross-section,? = S, spin, V given by: 

a px° = a ° (fs) 2 ( mjp ) 2 ( sca l ar ) > (the isovector scalar is negligible, i.e. 
°p = <£) 

= a 3 (/° + /A) 2 (^f) 2 (spin) , < x o = *o (/£ + #) a (^) 2 
(vector) 

where mjy is the nucleon mass, r\ = m x /mjv^4, and /x r (A) is the LSP-nucleus 
reduced mass, fi r (N) is the LSP-nucleon reduced mass and 

= ^-(G F m N ) 2 ~ 0.77 x KT 38 cm 2 (11) 

Z7T 



Q = Qo" , Qo = 1 , 9 = 4.1 x 10 4 A- 4 / 3 tfeF (12) 
Amjvo 

where Q is the energy transfer to the nucleus and F(u) is the nuclear form 
factor. 

In the present paper we will concentrate on the coherent mode. For a 
discussion of the spin contribution, expected to be important in the case of the 
light nuclei, has been reviewed elsewhere^. 



IV. Expressions for the Rates. 

The non-directional event rate is given by: 
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-K = Unon-dir = ~ TT = : <T{U,V)\v\ (13) 

at m x Aitln 

Where p(0) = 0.3GeV/cm 3 is the LSP density in our vicinity and m is the 
detector mass The differential non-directional rate can be written as 

dR = dR non ^ dlr = da(u,v)\v\ (14) 

m x ArriN 

where da(u,v) was given above. i— ir— i 

The directional differential rateliJpZl in the direction e is given by : 

p(0) m 1 

dRdir = : v.eHiv.e) — da(u,v) (15) 

m x ArriN Zir 

where H the Heaviside step function. The factor of l/2w is introduced, since 
the differential cross section of the last equation is the same with that entering 
the non-directional rate, i.e. after an integration over the azimuthal angle 
around the nuclear momentum has been performed. In other words, crudely 
speaking, l/(2ir) is the suppression factor we expect in the directional rate 
compared to the usual one. The precise suppression factor depends, of course, 
on the direction of observation. The mean value of the non-directional event 
rate of Eq. (|l4|), is obtained by convoluting the above expressions with the 
LSP velocity distribution J(v,ve) with respect to the Earth, i.e. is given by: 

idR\ p(0) m f da(u,v) 3 

The above expression can be more conveniently written as 

(f) = pM i^vw)(f) ,<£) = [ ^f(v,v E) d -^d 

\ du I m x ArriM du du J ./ (v 2 ) du 
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(17) 



After performing the needed integrations over the velocity distribution, 
to first order in the Earth's velocity, and over the energy transfer u the last 
expression takes the form 

R = Rt [1 + h(a, Q min )cosa)} (18) 

where a is the phase of the Earth (a = around June 2nd) and Q m in is the 
energy transfer cutoff imposed by the detectoirT-|In the above expressions R is 
the rate obtained in the conventional approachEj by neglecting the folding with 
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the LSP velocity and the momentum transfer dependence of the differential 
cross section, i.e. by 

P(0) 



R = ^^—V¥)[t s + t spm + ^ v ] (19) 

where Si, i — S, V, spin contain all the parameters of the SUSY models. The 
modulation is described by the parameter h . 

The total directional event rates can be obtained in a similar fashion by by 



integrating Eq. (15) with respect to the velocity as well as the energy transfer 
u. We find 

Rdir = R[(t di r/2ir) [1 + (hi - h 2 )cosa) + h 3 sina] (20) 

where the quantity t^ir provides the un modulated amplitude, while h\,h 2 
and /13 describe the modulation. They are functions of the angles Q and <E>, 
which specify the direction of observation e, as well as the parameters a and 
Qmin- The effect of folding with LSP velocity on the total rate is taken into 
account via the quantity tdir, which depends on the LSP mass. All other SUSY 
parameters have been absorbed in R. In the special case previously studied, 
i.e along the coordinate axes, we find that: a) in the direction of the sun's 
motion fi2 = /13 = Q, b) along the radial direction (y axes) h% = and c) in 
the vertical to the galaxy h 2 = 0. Instead of tdir itself it is more convenient to 
present the reduction factor of the un modulated directional rate compared to 
the usual non-directional one, i.e. 

fred = ^r = W(2?r t) - k/(2tt) (21) 

It turns out that the parameter k, being the ratio of two rates, is less dependent 
on these parameters. Given the functions hi(a, Q m in), ' = 1,2, 3, one can plot 



the the expression in Eqs (18) and E0 as a function of the phase of the earth 



V. The Scalar Contribution- The Role of the Heavy Quarks 



The coherent scattering can be mediated via the the neutral intermediate Higgs 
particles (h and H), which survive as physical particles. It can also be mediated 
via s-quarks, via the mixing of the isodoublet and isosinlet s-quarks of the same 
charge. In our model we find that the Higgs contribution becomes dominant 
and, as a matter of fact the heavy Higgs H is more important (the Higgs 
particle A couples in a pseudoscalar way, which does not lead to coherence). It 
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is well known that all quark flavors contribute Ha, since the relevant couplings 
are proportional to the quark masses. One encounters in the nucleon not only 
the usual sea quarks (uu, dd and ss) but the heavier quarks c.-k t which couple 
to the nucleon via two gluon exchange, see e.g. Drees et al and references 
therein. 

As a result one obtains an effective scalar Higgs-nucleon coupling by using 
effective quark masses as follows 

m u -> f u m Nl m d -> f d m N . m s -> f s m N 

rriQ — > /qtojv, (heavy quarks c,b,t) 

where m^r is the nucleon mass. The isovector contribution is now negligible. 
The parameters f q , q = u,d,s can be obtained by chiral symmetry breaking 
terms iixjrelation to phase shift and dispersion analysis. Following Cheng and 
Cheng E£l we obtain: 

/„ = 0.021, f d = 0.037, f s = 0.140 (model B) 

/« = 0.023, f d = 0.034, / s = 0.400 (model C) 

We see that in both models the s-quark is dominant. Then to leading order 
via quark loops and gluon exchange with the nucleon one finds: 

fq = 2/27(1 - E 9 /,) - i-e- fq = 0.060 (model B), f Q = 0.040 

(model C) 

There isr-a, correction to the above parameters coming from loons involving 
s-quarksl!2l and due to QCD effects. Thus for large tanft we findliil: 

f c = 0.060 x 1.068 = 0.064, f t = 0.060 x 2.048 = 0.123, f b = 0.060 x 1.174 = 

0.070 (model B) 

f c = 0.040 x 1.068 = 0.043, f t = 0.040 x 2.048 = 0.082, f b = 0.040 x 1.174 = 

0.047 (model B) 

For a more detailed discussion we refer the reader to RefsEHEl. 



VI. Results and Discussion 

The three basic ingredients of our calculation were the input SUSY parameters 
(see sect. 1), a quark model for the nucleon (see sect. 3) and the velocity 
distribution combined with the structure of the nuclei involved (see sect. 2). 
we will focus our attention on the coherent scattering and present results for the 
popular target 127 1. We have utilized two nucleon models indicated by B and 
C which take into account the presence of heavy quarks in the nucleon. We also 
considered energy cut offs imposed by the detector, by considering two typical 
cases Qrain — 0, 10 KeV. The thus obtained results for the un modulated non 
directional event rates Rt in the-ease of the symmetric isothermal model for a 
typical SUSY parameter choice t3 are shown in Fig. 
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140 




m LSP (GeV) 



FIG. 1.: The Total detection rate per (kg — target)yr vs the LSP mass in GeV for a typical 
solution in our parameter space in the case of 127 1 corresponding to model B (thick line) 
and Model C (fine line). For the definitions see text. 



The two relative parameters, i.e. the quantities t and h, are shown in Fig. 
|2| and Figs |^,^ respectively in the case of isothermal models. 




FIG. 2.: The dependence of the quantity t on the LSP mass for the symmetric case (A = 0) 
on the left as well as for the maximum axial asymmetry (A = 1) on the right in the case of the 
target 127 /. For orientation purposes two detection cutoff energies are exhibited, Q m i n = 
(thick solid line) and Q m i n = 10 keV (thin solid line). As expected t decreases as the cutoff 
energy and/or the LSP mass increase. We see that the asymmetry parameter A has little 
effect on the un modulated rate. 
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FIG. 3.: The same as in Fig. y for the modulation with A = 0. We see that the modulation 
is small and decreases with the LSP mass. It even changes sign for large LSP mass. The 
introduction of a cutoff Q m i„ increases the modulation (at the expense of the total number 
of counts) . 




FIG. 4.: The same as in Fig. |3j for A = 1. We see that the modulation increases with the 
asymmetry parameter A. 

i_. The case of non isothermal models, e.g. caustic rings, is more complicated 
E3 and it will not be further discussed here. 

It is instructive to examine tbe, reduction factors along the three axes, i.e 
along +z, —z, +y, —y, +x and — xc3. Since f re d is the ratio of two parameters, 
its dependence on Qmin and the LSP mass is mild. So we present results 
for Qmin = and give an average as a function of the LSP mass (see Table 
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. As expected the maximum rate is along the sun's direction of motion, i.c 
opposite to its velocity (—z) in the Gaussian distribution and +z in the case 
of caustic rings. In fact we find that «(— z) is around 0.5 (no asymmetry) 
and around 0.6 (maximum asymmetry, A = 1.0), It is not very different from 
the naively expected f ree i — 1/(2tt) = K = 1. The asymmetry \Rdir(~) ~ 
Rdir{+)\/ {Rdir{—)+Rdir{+)) is quite large in the isothermal model and smaller 
in caustic rings. The rate in the other directions is quite a bit smaller (see Table 



As we have seen the modulation can be described in terms of the param- 
eters hi, i = 1,2, 3 (see Eq. (^o|)). If the observation is done in the direction 
opposite to the sun's direction of motion the modulation amplitude h\ be- 
haves in the same way as the non directional one, namely h. It is instructive 
to consider directions of observation in the plane perpendicular to the sun's 
direction of motion (0 = 7r/2)even though the un modulated rate is reduced in 
this direction. Along the — y direction ($ = (3/2)ir) the modulation amplitude 
hi — h<x is constant, —0.20 and —0.30 for A = 0, 1 respectively. In other words 
it large and leads to a maximum in December. Along the +y direction the 
modulation is exhibited in Figs and ^. 




FIG. 5.: The quantity hi — h,2 in the direction +y for A = (thick line) and A = 1 thin line. 
In the —y direction this quantity is constant and negative, —0.20 and —0.30 for A = and 
1 respectively. As a result the modulation effect is opposite (minimum in June the 3nd). 
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FIG. 6.: The same as in Fig. KAfoT the modulation amplitude in the direction +x, which is 
essentially /13, since|hi| << \hs\. Thus in this case the maximum occurs around September 
the 3nd and the minimum 6 months later. The opposite is true in the — x direction. 



TABLE I.: The ratio k of the un modulated directional rate along the three directions to 
the non-directional one: z is in the direction of the sun's motion, x is in the radial direction 
and x is perpendicular to the axis of the galaxy. The asymmetry is also given. Q m i n = 
was assumed. 



isothermal I caustic rings 



A dir. + - asym + - asym 



z 0.02 0.50 0.92 0.75 0.25 0.50 

y 0.16 0.16 0.22 1.00 

x 0.16 0.16 0.37 0.24 0.21 

1 z 0.04 0.58 0.90 

1 y 0.12 0.12 - - - 

1 x 0.17 0.17 - - - 
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VII. Conclusions 



In the present paper we have discussed the parameters, which describe the 
event rates for direct detection of SUSY dark matter. Only in a small segment 
of the allowed parameter space the rates are above the present experimental 
goals. We thus looked for characteristic experimental signatures for back- 
ground reduction, i.e. a) Correlation of the event rates with the motion of the 
Earth (modulation effect) and b) the directional rates (their correlation both 
with the velocity of the sun and that of the Earth.) 

A typical graph for the total un modulated rate is shown Fig. |l|. The 
relative parameters t and h in the case of non directional experiments are 
exhibited in Fig. || and Figs || and [|. We must emphasize that these two 
graphs do not contain the entire dependence on the LSP mass. This is due 
to the fact that there is the extra factor /x^ in Eq. (2.10) but mainly due to 
the fact that the nucleon cross section depends on the LSP mass. Fig ^| and 
Figs H and || were obtained for the scalar interaction, but do not expect a 
very different behavior in the case of the spin contribution. The overall spin 
contribution, however, is going to be very different, but such considerations 
are beyond the goals of the present paper. We should also mention that in the 
non directional experiments the modulation 2h\ is small, .i.e. for A — less 
than 4% for Qmin = and 12% for Qmin = 10 KeV (at the expense of the 
total number of counts). For A = 1 there in no change for Q m in — 0, but it 
can go as high as 24% for Q m i n = 10 KeV. 

For the directional rates It is instructive to examine the reduction factors 
along the three axes, i.e along +z, —z, +y, —y, +x and —x. These depend 

a the nuclear parameters, the reduced mass, the energy cutoff Qmin an d A 
Since f re d is the ratio of two parameters, its dependence on Q m in and 
the LSP mass is mild. So we present results for Q m in — and give their 
average as a function of the LSP mass (see Table |). As expected the maximum 
rate is along the sun's direction of motion, i.e opposite to its velocity (—z) in 
the Gaussian distribution and +z in the case of caustic rings. In fact we 
find that k(— z) is around 0.5 (no asymmetry) and around 0.6 (maximum 
asymmetry, A = 1.0). It is not very different from the naively expected f rec i — 
l/(27r) = k — 1. The asymmetry along the sun's direction of motion, asym = 
\Rdir(—) — Rdir(+)\/ (Rdir(-) + Rdir(+)) is quite characteristic, i.e. quite large 
in the isothermal models and smaller in caustic rings. The rate in the other 
directions is quite a bit smaller (see Table ||). But, as we have seen, in the plane 
perpendicular in the sun's velocity we have very large modulation, which may 
overcompensate for the large reduction factor. It is interesting to note that 
the dependence on the phase of the Earth a changes substantially with the 
direction of observation. 

In conclusion in the case of directional un modulated rates we expect 
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characteristic and pretty much unambiguous correlation with the motion on 
the sun, which can be explored by the experimentalists. The reduction factor 
in the direction of the motion of the sun is approximately only 1 / (47r) relative 
to the non directional experiments. In the plane perpendicular to the motion 
of the sun we expect interesting modulation signals, but the reduction factor 
becomes worse. A more complete discussion will be given elsewhere t3. 

This work was supported by the European Union under the contracts RTN 
No HPRN-CT-2000-00148 and TMR No. ERBFMRX-CT96-0090. 
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